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Formulas are derived for the estimation of the variances of position lines from fixes 
with unknown target positions. Two approaches are considered, (1) that presented by 
Daniels [1951], and (2) an analysis of the squares of the errors in the position lines assuming 
the target is at the least squares estimate for its position. 

1. Introduction 

An important problem in position fixing is to estimate the accuracy of one's position 
lines. This can be done by assigning a variance to each line, see for example Daniels [1951] 
and Beale [1961]. These variances are usually estimated from fixes on targets whose true 
positions are known to the analyst. There are two major difficulties about estimating vari- 
ances from fixes where the true position is unknown: (1) The estimates depend critically on 
the assumption that errors in different position lines are independent; and (2) The statistical 
problem is difficult, and any valid method seems to involve a substantial computing effort. 

In spite of these difficulties, the problem deserves attention, because it is sometimes 
impossible to set up a satisfactory program of check fixes on known targets. Two approaches 
are considered in this paper. One follows Daniels [1951], the other is based on an analysis 
of the squares of the errors in the position lines assuming 1 hat the target is at the least squares 
estimate for its position. 

Section 2 of this paper presents the basic assumptions common to both these approaches. 
The next 3 sections are concerned with Professor Daniels' approach to. the problem. The 
approach is presented in general terms in section 3. Specific formulas for the important 
special case of 4 lines per fix are developed in section 4 ; and the application of these formulas 
is considered in section 5. The formulas for the alternative approach are developed in sec- 
tion 6, and their application is considered in section 7. 

Finally, in section 8 some suggestions are made for an artificial sampling experiment to 
try out both these approaches. 

2. Basic Assumptions 

The basic assumptions made in this work are: (1) The earth is flat near the true position; 
(2) the position lines are straight lines; (3) an error of observation displaces the line parallel 
to itself; (4) the errors have zero means and are statistically independent; and (5) the errors 
are normally distributed, and we have rough estimates of their variances. 

Given these assumptions, we can take Cartesian coordinates, and denote the j th position 
line by the equation 

x sin Bj—y cos B j =p il (2.1) 

where 6j is a known constant, and pj is regarded as a random variable with mean 

£ sin 6j— rj cos dj, 



1 Contribution from Admiralty Research Laboratory, Teddington, Middlesex, England. 

2 Paper presented at the Conference on Transmission Problems Related to High-Frequency Direction Finding, at UCLA, June 21-24, 1960. 
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and variance a), where (£,77) are the unknown coordinates of the true position. In our theo- 
retical work, we assume that £=77=0, so that the mean value of pj is zero; but we must re- 
member that pj is then not a directly observable quantity. 

Some discussion of these basic assumptions may be helpful. 

The first 3 assumptions are introduced to make the model linear, i.e., to make the value 
of each observation in the absence of errors a linear function of the unknown parameters 
£ and r). It is clear that in some DF situations, notably when the DF stations have a narrow 
base line, the problem is decidedly nonlinear when expressed in terms of obvious parameters, 
such as the distances of the target east and north of some origin. But there is evidence to 
suggest that the "intrinsic nonlinearity" as defined by Beale [1960], is nearly always small 
for DF problems. This implies that any results derived from the standard linear theory will 
be valid, provided that they are expressed in terms independent of any particular system 
of parameters (and provided of course that any other assumptions made about the error 
distributions are valid). 

Using the language of Beale [1960], we must define all the quantities we use in terms of 
sample space and the solution locus, and not solely in terms of parameter space (or possible 
target positions on the Earth, which is a representation of parameter space). Now a position 
line can be defined in sample space as the intersection with the solution locus of the hyperplane 
where one coordinate, i.e., one observation, is constant. All we require in this work is an inter- 
pretation of the quantities p jf jf and a). This can be obtained by taking some specific point T 
on the Earth near the true target position, and imagining a parameter system that coincides 
locally with distances east and north of T. Now let /3j denote the actual ,7 th observation. (In 
the DF problem this is the bearing from the j th DF station.) Let f3 jT denote a hypothetical j th 
observation giving a position line passing through T, and let pj(ffj) denote the (signed) distance 
on the Earth of the j th position line from T as a function of (3j. 

Let \j denote d Pj/dfij, evaluated where Pj=/3 jT . 

Then we can write pj=\ j(Pj-Pjt), &)— X^ var (Pi)> an d ^ * s the angle that the hypothetical 
,7 th position line passing through T makes with some arbitrary ar-axis passing through T. 

Note that it is not necessarily legitimate to interpret Bj as the angle that the observed 
position line on the Earth makes with some 2-axis; or to interpret Xj as dpj/dftj evaluated for the 
observed fij. It is important that the same point T should be used for all position lines in the 
fix. (This approach can be used to derive the standard (Gauss) iterative procedure for finding 
the least squares estimate for the target position when the variances of the (3j are assumed 
known. The point T must then be taken as the trial value of the least squares point at each 
stage.) 

Our 4th basic assumption, that the error distributions have zero means and are inde- 
pendent, is very important. In certain circumstances it will not be satisfied in practice, in 
which case the least squares point will usually be a less accurate estimate of the target's true 
position than standard theory suggests, even if the variances are correctly estimated; and 
furthermore the variances will be underestimated by the methods described in this paper. If 
we had numerical values for all correlations involved, we could allow for them, but in general 
there seems to be no practical alternative to assuming independence, hoping for the best, but 
realizing the possibility of being misled by correlated data. 

Our 5th assumption, that the errors are normally distributed and we have rough estimates 
of their variances, is not so critical. It is only required to give appropriate weights in our least 
squares estimation of the variances, and to estimate the accuracy of our final variance estimates. 

3. Daniels' Approach in General 

The problem considered in this paper was considered by Daniels [1951]. In section 9, 
Daniels writes: 
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"When the true position corresponding to each cocked hat is unknown, ... it is found that the method of 
maximum likelihood is useless for estimating the variances of the lines. It gives in general inconsistent answers, 
and sometimes no answer at all ... . The method [breaks down because] the number of incidental parame- 
ters to be estimated increases with the number of cocked hats used. 

"Nevertheless a simple method of estimating the variances does exist requiring no assumptions about the 
error distributions, provided parallel displacement is assumed and the errors are independent. Suppose there 
are 2V three-line cocked hats, not necessarily for the same true position, the lines having variances <rf, a% a\ 
to be determined. The expression 

u = pi sin (0 2 — 3 )4-p2 sin (0 3 -0i) + Pz sin (0i-0 2 ) (3.1) 

has the same value whatever the point from which the perpendiculars p h p 2} and p3 are dropped on to the lines, 
and u is easily measured by taking the point at a vertex. ... If the point is chosen at the unknown true po- 
sition it follows that 

E(u 2 )=*l sin 2 (02-0 3 )+<ri sin2 fa-ej+al sin 2 (0i-0 2 ). (3.2) 

If iV>3, the N equations with the observed values of u 2 replacing E(u 2 ) may be solved by least squares. Since 
for normal errors var (u 2 ) = 2{E(u 2 ) } 2 , arbitrary weights have to be used in the first solution and correct weights 
approximated to in subsequent iterations. The estimates are, however, unbiased for any error distribution 
whatever the weights used." 

But it should be noted that with fixed I)F stations the normal equations derived from (3.2) 
to estimate the variances are liable to be very ill conditioned, and may even be singular. If the 
simplifying* assumption about parallel displacements were strictly correct, then the values of 
0i, 2 , and 3 would be the same for all tasks, so that the right hand sides of all equations of the 
form (3.2) would be strictly proportional, and one could never hope to estimate more than the 
given linear function of erf, <j 2 , and of. One 4 might hope in practice to be saved by the error in 
the approximation; since d 3 may be effectively constant for all reasonably possible observa- 
tions on a given target, but not constant over all targets. But this will only be so if the targets 
are well distributed about the DF stations. 

It is fairly easy to see that this difficulty is fundamental to the problem, and is not simply a 
defect in the present approach. For if the 0, were strictly constant, one could for example 
assume that o- 1 = o- 2 =0, and the signed distance of the point of intersection of the first two po- 
sition lines from the third will have a certain probability distribution (normal if the bearing errors 
are normal) with mean zero and variance given by (3.2). This variance is constant as long as 
the true values of a'i, ai, and o\ are constant. So data of this type can never disprove the 
hypothesis that ai = a 2 = 0. 

In the analogous one-dimensional problem, if instruments 1 and 2 give independent obser- 
vations Xi and x 2 of an unknown scalar quantity £, and x t has mean £ and variance tr?, then 
E(x x — x 2 ) 2 = g\-\- o-\. With an arbitrarily large number of pairs of observations one can therefore 
get an arbitrarily accurate estimate of o-i + o- 2 , but cannot estimate a 2 or <j\ separately. But if a 
third independent instrument is available, then one can estimate <r 2J rcr 2 2 by (xi— x 2 ) 2 , a\ 
+ 0"! by (xi—Xs) 2 , and ci+ci by (x 2 — # 3 ) 2 . Hence a 2 is estimated by the mean value of h(xi—x 2 ) 2 
+ § ( x i — ^3) 2 — i (^2 — ^3) 2 , i.e., x\— X\X 2 — xix 3 +x 2 x 2 , as pointed out by Pearson [1902]. 

If one adds a fourth position line to the 2-dimensional problem, then one can obtain four 
expressions analogous to (3.2) by taking each set of three lines in turn. It turns out that these 
four expressions do not suffice to estimate the four variances o-f, a 2 , 03, and a 2 ; because the deter- 
minant of the left hand side of the four equations vanishes identically. But Professor Daniels 
has pointed out in a private communication that one can get more expressions by considering 
tin 4 mean values of the products of ^-statistics for different triangles. Thus if 

u m = pi sin (0 2 -0 3 ) + p 2 sin (0 3 — 0i) + p 3 sin (d l — 2 ), 
and 

U 12i =p l sin (0 2 -0 4 ) + p 2 sin (0 4 — Oj+pt sin (0i~0 2 ), 
then 

E(u 123 u V2i ) = (T l 2 sin (0 2 — 3 ) sin (0 2 — d A ) + a 2 2 sin (0 3 — 0i) sin (0 4 — 0i). 
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Figure 1. A possible set of 4 position lines. 

Attention must be paid to signs. It is convenient to regard each position line as a directed 
line, pointing in the direction 6 from some apex at an arbitrarily large distance from the scene 
of action. One can then define the corresponding displacement p as positive for points whose 
bearing from the apex is greater than 6 (by an arbitrarily small amount) and negative on the 
other side of the position line. This will ensure that the signs of the terms on the right hand 
sides of the expressions analogous to (3.1) come out as expected. But it will still be necessary 
to study the configuration of lines to see which ^-statistics are positive and which are negative. 
In the situation illustrated in figure 1, u iu and ^34 are positive, and ^123 and u 12 ± are negative. 

Given fixes with n^>3 position lines per fix, one can obtain U=n(n-l)(n-2)/6 u-statistics 
from each fix, and U(U-\-l)/2 derived statistics by considering the squares and products of 
the w-statistics. 

If n =4, this gives 4 ^-statistics and 10 derived statistics. 

If n=5, we have 10 u-statistics and 55 derived statistics. 

This approach is therefore not very practical in its present form for n>5. But the case 
n=i is important, as it is the smallest n for which unique variance estimates can be obtained 
from a set of fixes, each with essentially the same values of 6 j} though not all on the same tar- 
get. We therefore explore this case in more detail in the next section. 

4. Daniels' Approach With n==4 

With 4 position lines one can form 4 triangles, and obtain 4 basic quantities u m , u n ^ u lu , 
and u 2 34, which have the following expressions 
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^123 = 01*23 + 02*31 + 03*12 

^124 = 01*24 + 02*41 + 04*12 

^134 = 01#34 + 03#41 + 04*13 

^234 = 02#34 + 03#42 + 04#23, j 

where s tj denotes sin (fit— &j). 

From these one can derive a vector y of 10 derived observations, such 

that 
where 



and A is 
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where ti=s l2 s M , t 2 = s 13 s 2i and t 3 =s H s 23 . 
Then 



E{y) =-&(*?, d, 4, <%', 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



where II is the matrix formed from the 1st, 5th, 8th, and 10th columns of A, since the pj are 

independent and have zero means when referred to the true target position as origin. 

But in order to determine rational, as opposed to arbitrary, least squares estimates for the 
unknown parameters <r\ } <j\, 0-3, and o-f, it is necessary to consider the co variance matrix of y, 
We rely on the fact that the pj are independent and have zero means, and find that 



where 



c iJ =cov{y i y j ) =Ey i y j — (Ey t ) (4)=S«<W, 

k 

2 4 9 9 9 9 9 9 "\ 

01, r 2 =at(T2, r 3 =<rf(7i, r 4 =o-io-4 " 

r 5 =2at, r & =a 2 2 al r 2 =<r 2 2 <rl 

r 8 =2at, r 9 =al<rl 

r 10 =2a\. 



r 



(4.7) 



(4.8) 



(The coefficients 2 in r h r s , r 8 , and r 10 should be replaced by 2 + k 4 if the error distributions have 
a fourth cumulant k 4 different from zero.) 
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The co variance matrix C=(c i3 ) of the derived observations is therefore obtained from 
(4.7), or equivalently by defining a new matrix 

where f ij =a ij r) /2 , 

sothat C = FF'. 

Now the standard least squares results can be expressed in matrix notation as follows : 
If Y=X/3+e, where Y is a vector of observations, X is a known matrix, f3 is a vector of unknown 
parameters, and e is a random vector whose components are uncorrelated and have means 
zero and a common variance o- 2 , then the least squares estimates B of /? are given by choosing 
B to minimize 

Q=(XB-Y) , (XB-Y), (4.9) 

Assuming that the errors are normally distributed, the likelihood is proportional to 
exp (—Q/2a 2 ), and in any case 

B=(X , X)- 1 X / Y, (4.10) 

and the co variance matrix for the estimates B is given by (X' 'K)~ 1 a 2 . (4.11) 

If now y=GY, we have y=GX/3+Ge, and therefore the co variance of y t and yj is given by 

k 

Hence we can identify G with the matrix F, putting a 2 =l, and we know that, if /3=(a 2 , 

(7 2 , 0% V±) , 

y=H£+Fe. 

So X=F- 1 H, J=F"V, 

and 

Q={F-i(H£-£)}'{F-'(H£-y)} 

= (Hg-y)'F" 1 'F- I (H^-y) 
= (H^-y)'(FF')-'(Hg-y) 

= (H£-y)'C-(H£-y). (4.12) 

Further 

£= (H / F / - 1 F- 1 H)- 1 H / F / - 1 F" 1 y 

= (H / C- 1 H)- 1 H / C- 1 y, (4.13) 

and the covariance matrix of is given by 

(H'C" 1 !!)- 1 . (4.14) 

These results are due to Aitken [1935]. 

5. Application of the Formulas of Section 4 

To apply these formulas to a situation in which one has N fixes, i.e., sets of four position 
lines, each with a constant variance and a constant value of 6 h one proceeds as follows. 
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One first sorts out one's sign convention, as indicated at the end of section 3. This will 
produce signed values of u 12 z, u^, u m , and u 2 s4 for each of the Affixes. From these one computes 
derived observations y h . . ., ?/ 10 for each fix from (4.3). Averaging over all fixes, one obtains 
a vector 

y=(»i, $2, • • ■, 2/io)'. 

We next compute the matrix A from (4.5). 

By using the guessed values for the unknown variances cf, er|, erf, and erj, one can derive 
values of r 2 , r 2 , . . ., rf from (4.8), and hence compute the elements r,, of (lie matrix C 
from (4.7). 

This matrix must be inverted to produce C _1 . 

We must also consider the matrix II whose columns are the 1st, 51 h, 8< h, and 10th columns 
of A, and our estimates of a'i, ai, o-§, and <r\ are the 4 elements of the matrix 

(H'C^H^H'C" 1 ^ (5.1) 

If our guessed values of the of, used to define r, were correct, the covariance matrix for 
these estimates would be 

l(H'C-'H)-'. (5.2) 

One could try iterating this procedure, using the estimated variances to produce a new 
vector r, and hence a new C, and hence revised estimated variances, and so on. 

1 am doubtful of the wisdom of this, though it might be interesting to try it on an artificial 
sampling experiment. If one wants a reasonable estimate of the accuracy of the estimates one 
must compute a revised ( 1 to use in (5.2). But I can see little merit and some danger in revising 
the weights to produce revised estimates: Note that the iterated estimates are not unbiased. 

This procedure may on occasion produce negative estimated variances. Steps must be 
taken at least to ensure thai negative (or even wry small positive) variances are not used in 
the definition of r. Theoretically the least squares problem should be set up as a quadratic 
programing problem when such negative estimates are produced by the standard method. 
But this is too large an issue to face at this stage. 

6. A Direct Approach 

In this section we derive some equations concerning the apparent error in the position 
lines, assuming that the target is at the least squares estimate for its position. 

Suppose that the j th position line is 

x sin Bj—y cos dj=Pj, 

with the origin at the (unknown) true target position. Suppose further that we have some 
rough estimate s 2 for a 2 . 
We write 

(7 ? = ^y(/l+^) 7 an d l/sj=Wj. 

Now r the sum of squares function, whose minimum value defines the least squares 
estimate, is 



where 



Q=^2wj(x sin Bj—y cos Bj—pj) 2 =a x 2 -\-2hxy J rby 2J r2gx J r 2fy J rc, 



a=^P,w j sin 2 B h b = ^2wj cos 2 B h c — y^p 2 . 



j 



9 } .) 



f=^LjWjPj cos B h g= —^,Wjpj sin h h= — S^j sm Oj cos 0j> { 
3 j 
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Note that a, b, and h are known constants, while /, g, and c depend on the random 
variables pj. 

It is convenient to rotate axes, so that in the new coordinates (X, Y), 

Q=AX 2 +BY 2 + 2GX+2FY+C. 

This involves finding an angle a such that 

fa=dj+a, (6.2) 

X sin fa—Y cos fa=pj, 

X sin (dj+a)— Y cos (Oj+a)=pj, 

x=X cos a-\-Y sin a, 

y=—X sin a+F cos a. 

And since ax 2 +2to/+ 6?/+2^+2y7/+c=JLZ 2 +2i7Zr+BF 2 +2G ? A : '+2^r+(7 J where H=0, 
we have 

^l = a cos 2 a—2h cos a sin a+b sin 2 a, 

0=a cos a sin a-f A(cos 2 a— sin 2 a) — b cos a sin a, 

B=a sin 2 a+2A cos a sin a-\-b cos 2 ck, 

#=# cos a— /sin a, 

i^=0 sin a+/ cos a, 

(7=c. 



and the jtii position line is 

i.e., 

so that 



(6.3) 



This implies that 



tan 2a=2h/(b — a). 



(6.4) 



This defines a set of values of a differing by multiples of 90°. It does not matter which 
is taken. 

The values of A, B, C, F, and G can be obtained either from (6.1) and (6.3), or alter- 
natively from the formulas 



(6.5) 



A=^Wj sin 2 fa, B=^2wj cos 2 fa, C=y^,w 4 ph 

j j 3 

F= y^.Wip, cos fa, G=—^2w j p j sin <j>j, H=—^2wj sin fa cos <^=0 ; 

j j j 

where <t>j is obtained by substituting for a from (6.4) in (6.2). 

The transformed coordinates of the least squares estimate are then 

X =-G/A and Y =-F/B. (6.6) 

Now let d k denote the signed displacement of the k th position line from the least squares 
estimate. This quantity can be observed in practice. Then 



where 



G F n 

d k =Pk+^ sin fa— -g cos fa=^ PiMw k ) 1/2 \ ik , 



t—w k 



sin 2 fa c os 2 fa 
+ B 



> 



Note that 



x . , m/o /sin<fo sin fa cos fa cos fa\ ,., 7N 

h k =hd=-(WiWk) m y ^ ~-\ ^g -)> {%^k). 



S \ub=ft— Sw* sin2 W^-S^* cos2 falB=n—2. 

k=l 



(6.7) 

(6.8) 
(6.9) 

(6.10) 
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Further 

COS <j>i Sin <t>k C0S <l>k 



» 2 _ / sin 2 fa , cos 2 ^ \ , » f w, sin 2 0,- sin 2 <fa , w f sin <fo 



AB 



So 



-0^ cos 2 <fo\ - /sin 2 <fe cos 2 </>A 

"5" "/"'"^V-T" 4 --^-/ 



S X?*=X». (6.11) 



Further, if ^f^Z, 

X^\ ^ o\ i /„.. ... m/2 x^ ... / sin & sin< ^ i cos </>* cos #A / sin ^ sin ^ I (>QS ^ CQS M 

ox i / N1/2 /sincfe sin</>, cos 0* cos </>A „ , N 
=2X^+(w*<) 1/2 ( — ^ — ~H — ^g— : ) from (6.5) =X W . 

So (6.11) can be generalized to read 

£xoX<i=\ki. (6.12) 

*=1 

From (6.7) we deduce that 

w*rf*=S S {WiWjY'^ik^jkPiVj- (6.13) 

i = l j'=l 

Now 

#2>?=ff?=(l+S 4 )M, Ep tPi =0 (i*j), 

So 

Sw*di=S X?t(l+5,)=X«t+S «A?», from (6.11). (6.14) 

i=l i = l 

It follows that 

E 2 w k dl=Tj X**+2 «/ 2 X?*=n-2+2 8i\a, from (6.10) and (6.11). (6.15) 

k=l k=l i=l &=1 i=l 

If 5^=5 for all i, we have the well-known formula 

ff&fcdi=(n-2)(l+*), (6.16) 

(which can easily be deduced from first principles). 

It is also of interest to consider the covariance of w k d\ and w t df. 
We have 

n 

cov (w/£, Wjd?)=S w?X?*X?,2(l+5,)7w? 

+222 v) i w j \ i1 ^ jl ^ i {k i i(X+b^{l + b j )lw i v) j by analogy with (4.7) 

i = l j^i 
n n 

=2 2 2 X rt X it X iI X i/ (l + 5 i )(l+«y). (6.17) 

i = l j=l 

It is of particular interest to consider the covariance matrix when all 8i=0, since this 
corresponds to the originally estimated covariance matrix. 
Using (6.12) we find this reduces to 

cov (w/£,wi?) = 2XL. (6.18) 
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If n=3, it can be shown that Xii/A 2 i=Ai2A22=Xi3A23> etc - ^ or example, to prove that 
AiiA22 — A? 2 = 0, one substitutes for the \ tj from (6.8) and (6.9) and multiplies out. The terms 
in 11 A 2 and IjB 2 cancel, so the expression can be written as a quadratic function of w u w 2 and 
w z (using the definitions of A and B) divided by AB, and II proves to be a factor of the 
numerator. 

This is an algebraic proof of something we already know on statistical grounds: that with 
only 3 bearings the random variables Widf, w 2 d\ and w z dl are totally correlated, and it is 
impossible to estimate more than one linear function of the unknown variances. 

7. Application of the Formulas of Section 6 

To apply the formulas of section 6 to a situation in which one has Affixes, i.e., sets of n 
position lines, each with a constant variance and a constant value of Q h one proceeds as follows. 

One starts with some assumed variances s) for the position lines, and one uses those to 
compute a least squares point estimate for the target for each fix. One then computes Wjd), 
i.e., d 2 /s 2 , being the weighted square of the displacement of the j th p ositio n line from each fix. 

One then averages over all N fixes, obtaining n quantities Wid\, . . ., w n d 2 n . 

One must also compute \ u for i=l, . . ., n, and \\j for i, j=l, . . ., n. These quantities 
are the same for all fixes. 

Note that \ it is defined by (6.8) and \ tJ by (6.9) if i^j. 

It may be more convenient to use the formula 

A2 _ (^^sin 2 z -)(^sin 2 i ) , 2(w i sm<focos0 < )(wjsin0jcos0 j ) , (w t cos 2 <fa) (wj cos 2 fy) r ~ ( . 
Then, from (6.14), we solve the equations 



S»«Mi=t«*3-\*, (t=i,...,»). (7.2) 

2 = 1 

and the resulting values of 8 t define the estimated variances, since 

<r«*«*« 2 (l + $«). 

To find the estimated co variance matrix for these estimates 8 U we form the matrix H 
such that 

and the matrix C such that 

(c ij ) = 2\ ij 2 (from (6.18)). (7.3) 

The estimated co variance matrix is then given by (5.2), if the 8 t prove to be small, and 

4 
this reduces to ^ C l . 

If one wants to combine these data with other data, then the contribution of these data to 
the sum of s quare s to be minimized is then NQ., where Q is given by (4.12) with 3=(8 l , . . ., 
8 n ) f and y=(w 1 di 2 , . . . ,w n d n 2 )'. 

8. Proposed Computational Program 

Some computational experience with both the schemes proposed here would be very 
valuable. 

In the first instance, it would be of interest to inspect the co variance matrix (4.14) for 
estimates obtained by Daniels' approach, i.e., with H defined following (4.6) and C defined 
by (4.7), for various values of 6j and a 2 . 

Ill 



1 1 would also be of interest to inspect the corresponding matrices for the direct approach, 
i.e., with (hij) = (X,/) and (c tj ) = (2X</). In this case (4.14) reduces to 4 C" 1 . The (ij) th element 
of this co variance matrix for the 8 must be course of multiplied by s*s* to represent the co vari- 
ance matrix for the variances. It will then be directly comparable with the other. 

One may hope to throw considerable further light on these procedures by an artificial 
sampling experiment. One must first choose values of 6 j} <r*, and (for the direct approach) 
8*. These will define values of a, b, and h from (6.1), and hence of A, B and a from (6.3) and 
(6.4). One can then compute the \ ik from (6.9). 

One then takes sets of values of p t as pseudo-random normal deviates with means zero 
and variances c* 2 , and use these to generate synthetic ^-statistics from (4.1), and values of d t 
from (6.7). One can then compare the estimated variances obtained by the two approaches 
(a) after one iteration, and (b) after using the results of each iteration as starting values for 
the next (i.e., to compute the r* for Daniels' approach and the w t for the direct approach) 
and continuing until the variances that go in come out. 

One can also consider the effects on the two approaches of using poor starting approxi- 
mations to the variances. 

It seems likely on general grounds that Daniels' approach should give better results, since 
it uses more information. But if the direct approach is at all satisfactory with n = 4, one may 
hope that it will he still better with larger values of n. 
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